Behavior of Ising Spin Glasses in a Magnetic Field 
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We study the existence of a spin-glass phase in a field using Monte Carlo simulations performed 
along a nontrivial path in the field-temperature plane that must cross any putative de Almeida- 
Thouless instability line. The method is first tested on the Ising spin glass on a Bethe lattice 
where the instability line separating the spin glass from the paramagnetic state is also computed 
analytically. While the instability line is reproduced by our simulations on the mean-field Bethe 
lattice, no such instability line can be found numerically for the short-range three-dimensional model. 
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Since its proposal in the mid-70 's, the Edwards and 
Anderson Ising spin-glass Hamiltonian has become a 
source of inspiration in statistical physics, especially in 
the context of mean field theory [l], y, 0] , and has been 
applied to a wide variety of problems across scientific 
disciplines. However, basic yet simple questions about 
the very nature of the spin-glass state in (experimentally 
relevant) finite space dimensions are still subject of con- 
troversy. The most prominent of such open questions is 
the existence of spin-glass ordering in a magnetic field. 

The fully-connected mean-field version of the EA 
model, called the Shcrrington-Kirkpatrick (SK) model 
, was solved using the replica method by Parisi [!, 0| ; 
the obtained free energy recently proven to be rigor- 
ously exact by Talagrand In this model, the low- 
temperature spin-glass phase is characterized by a com- 
plex free energy landscape made of many different val- 
leys. This phase also exists for low externally-applied 
magnetic fields and the so-called de Almeida-Thoulcss 
(AT) instability line [|| separates the spin-glass from the 
paramagnetic phase at finite ficlds/tempcratures. A sim- 
ilar scenario arises in the Bethe lattice approximation 
(see Fig. [I]). However, within the more phenomcnological 
description of spin glasses known as the droplet picture 
[?[ any infinitesimal field destroys the spin-glass order, in 
stark disagreement with the aforementioned mean-field 
description. Numerical evidence favoring the absence of 
a spin-glass state in a field for short-range spin glasses 
below th e upper critical dimension have become stronger 
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Capitalizing on the success of studying the suscepti 



bility and the finite-size correlation length [2CJ, l2l( to 
probe the spin-glass phase, we study the problem using 
a novel numerical approach backed up with analytic cal- 
culations on the Bethe lattice. On the numerical side, 
we use a multi-spin coded version of exchange (parallel 
tempering) Monte Carlo [13] with a new twist where the 
replicas "live" in the H-T plane along a nontrivial path 
that guarantees a crossing with a potential AT line (see 




FIG. 1: (Color online) Magnetic field H versus temperature 
T phase diagram of the Ising spin glass on a mean-field ran- 
dom regular graph with connectivity c = 6 with Gaussian 
interactions computed using the cavity method. The para- 
magnetic (PM) phase is separated from the spin-glass (SG) 
phase by the de Almeida-Thouless line (light red curve). The 
diagonal (dark blue) line h(T) represents the simulation path 
followed in the Monte Carlo simulations. Because the path 
starts at T > T c for H = and then increases along a diag- 
onal h(T) = r(T - T min ) in the H-T plane above T min < T c 
up to a value beyond T c (dashed vertical line) , an intersection 
with a putative AT line is guaranteed. The same approach is 
used for the diluted 3D simulations where T c — 0.663(6) [l9l ]. 



Fig. [1} and that in contrast to Refs. [23j, l2J] has both 
ends in the high temperature region where decorrelation 
is fast. We solve the thermalization problems on diago- 
nal paths found in Ref. 1(| by studying the link-diluted 
version of the model [19 [where the cluster moves intro- 
duced in Refs. [25| and 19[ can be used. In this model, 
the cluster updates allow the thermalization times to be 
decreased by a factor of at least 10 3 , allowing us to probe 
relatively large system sizes down to very low tempera- 
tures and large magnetic fields. We first demonstrate the 
efficiency of our strategy by applying it on the model de- 
fined on a regular random graph — which corresponds to a 
Bethe lattice in this context — for which we compute the 
AT line analytically, generalizing the original result by de 
Almeida and Thouless for the SK model Q . The analyt- 
ical and numerical results for the Bethe lattice agree to 
high precision. However, results on the 3D model show 



2 



no sign of an AT line. The fact that we do observe an 
AT line for the mean-field Bethe lattice, exactly where 
it is predicted, and do not for the short-range 3D Ising 
spin glass is convincing evidence that the phase diagram 
is indeed trivial and there is no spin-glass state in a field 
in three dimensions (3D). 

Models - The spin-glass Hamiltonian is given by 



H 



Jij &i Sj H Si . 



(1) 



The Ising spins Si £ {±1} have nearest-neighbor interac- 
tions. We study the mean- field case where N spins lie on 
the vertices of a regular random graph, and the EA model 
on a cubic lattice of size N = L 3 with periodic boundary 
conditions, where the interactions Jij £ { — 1,0,1} are 
chosen from a link-diluted bimodal distribution with a 
link occupation probability of 45% (l9j . 

Theoretical predictions - Following recent progress 
in the study of finite-connectivity mean-field systems, it 
is now possible to study (quite) precisely spin-glass mod- 
els on Bethe lattices using the cavity method [26j . Within 
this formalism, first steps in computing the phase dia- 
gram in the field-temperature plane have been achieved 
(see Refs. [27( and [14|). In order to determine the AT 
line, one needs to compute the onset of divergence of the 
spin-glass susceptibility \ in the paramagnetic phase in 
which case, thanks to the finite correlation length, the 
model is equivalent to a spin glass on an infinite tree 
with random boundary conditions [2(| (this is reminis- 
cent of the exact analysis presented in Ref. [2a|). We 
thus consider here the simplest Bcthc-Pcicrls cavity ap- 
proach for the problem and compute the AT line in this 
context. We refer the reader to Refs. (26[, 27 1, [T3|, and 
[29l | for a description of the method. For H = 0, we re- 
cover the well-known result first found by Thouless [30I ]: 
fc[tanh(/? c J) 2 ] av = 1, where (3 C = 1/T C , k = c - 1 (c the 
connectivity), and [■ • ■ ] av is a disorder average. For finite 
external field, one has to perform a numerical solution of 
the cavity equations to solve the cavity recursion and to 
compute the point where the susceptibility diverges [29| . 

Using this approach, we compute points along the AT 
line for different values of the field H with Gaussian 
and bimodal disorder and connectivities c = 3 and 6. 
We have also checked that the large-connectivity limit 
of our computation yields the SK result. We find that, 
close to T c , the data scale as h AT {T) ~ (T c - 7 1 ) 3 / 2 — 
as first predicted in Ref. and checked in Ref. [271]— 
and the instability line is then approximately linear close 
to zero temperature [HI]. A fit to /iat(7 1 ) = a\{T c — 
T) 3 / 2 exp(a2T + a 3 T 2 ) gives very accurate results whose 
precision is within the error of our numerical evalua- 
tion of the cavity recursion. For c = 6, we obtain 
01 = 0.786 (0.875), a 2 = 0.111 (0.221), and a 3 = -0.054 
(-0.127) and T c = 1.807 (2.078) for Gaussian (bimodal) 
distributed disorder. For c = 3, we find a\ = 0.785 
(0.827), a 2 = 0.251 (0.413), and a 3 = 0.117 (0.083) and 



T c = 0.748 (1.135). The Gaussian case with c = 6 is 
shown in Fig. [1] (light red curve). 

The presence of the instability line in finite- 
dimensional systems has been criticized before, especially 
in the context of the droplet model 0]. Other pictures 
were proposed where no such line is present [32I l33l. [34| . 
It has been also suggested that the line disappears below 
the upper critical dimension d u = 6 El, EI Sill, and 
even more complex scenarios are possible [3S|, |39|, |40| . We 
now attempt to clarify this issue in the 3D case. 

Numerical method — The simulations are performed 
using exchange Monte Carlo (EMC) (22)]. Traditionally, 
the field is fixed at a constant value and replicas at differ- 
ent temperatures perform a Markov chain in temperature 
space. There have been alternate approaches where the 
temperature T is fixed to T < T c and the replicas per- 
form a Markov chain in field space. This method has not 
proven to be efficient because tunneling across H = 
requires special moves [23j, [24| and especially because in 
this case no replica has the chance to reach the zero-field 
high-T phase where relaxation is fast. In order to solve 
this problem, we propose an EMC method in the H— 
T plane (see Fig. [1]) at which part of the replicas have 
labels (T, H = 0) for temperatures values in the range 
Tinin < T < T max with r max » T c (in 3D T min = 0.30). 
We also couple the replicas in the EMC scheme to a sec- 
ond set of replicas along a diagonal path in the H-T 
plane with labels (T, H) and h(T) = r(T - T min ), r con- 
stant. It is important not to choose the slope r of h(T) 
too steepbecause EMC is least efficient for large fields 
[H, [H, HH . In addition, T max (H > 0) > T C (H = 0) 



has to be chosen to ensure that the a potential AT line is 
crossed. The method has the advantage to deliver data at 
zero (horizontal blue path in Fig. [1]) , as well as finite field 
(diagonal blue path in Fig. [1]) . Finally, it is advantageous 
to choose the slope of the path in the H-T plane such 
that it crosses the phase boundary orthogonally (cleaner 
signal of the transition). 

We test equilibration by a logarithmic binning of the 
data until the data in the last three bins agree within 
error bars. To access relatively large system sizes in 
3D, we use a multi-spin coded version of the program 
that updates 32 copies of the system in parallel. In or- 
der to obtain thermalization in reasonable time, it is 
necessary to apply the cluster algorithm introduced in 
Refs. [19( and [25j. The speedup obtained by this ap- 
proach is of three orders of magnitude over conventional 
approaches, even for the smallest sizes simulated. Note 
that the speedup increases with increasing system size; 
our simulations would not have been possible otherwise. 
Simulation parameters are listed in Table [U 

Results on the Bethe lattice We have com- 
puted the connected spin-glass susceptibility using \ = 
N [(<1 2 )t - {q) 2 T ] av where q = l/Nj^Zi^ with U = 
SjSf being the overlap at site i between two indepen- 
dent replicas of the system. Here (• ■ ■ )t represents 
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TABLE I: Simulation parameters: iV sa is the number of sam- 
ples, iV sw is the number of Monte Carlo sweeps for one sample, 
T m i n is the lowest temperature simulated, -ff max is the max- 
imum field studied and N T is the number of replicas used in 
the EMC method. 7ibl(T) = 1.0(T-0.5) for the Bethe lattice 
(top set ); in 3D h w (T) = 0.7(T - 0.3) (bottom set). 
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FIG. 2: (Color online) Rescaled susceptibility x/N 1 ' 3 from 
Monte Carlo data of the regular random graph model with 
Gaussian disorder and connectivity c = 6. Left: Data for zero 
field cross at the analytically predicted critical temperature 
T c = 1.807. Right: Data along hsh(T) = 1.0(T - 0.5). The 
crossing point agrees very well with the values in Fig. [T] 



a thermal average. The susceptibility of the mean- 
field model has a finite-size scaling form x(T, H) = 
Ari/3<g (iv 1 / 3 ^ - T C (H)}) [H, [H hence data for differ- 
ent N should cross at T c when plotted as x/N 1 / 3 versus 
T. This is shown in Fig. [3] for c = 6 and a Gaussian 
distribution of the interactions: The left panel is at zero 
field, whereas the right panel is along the diagonal path 
^bl(T) = 1.0 (T — 0.5). In both cases, the data cross 
in agreement with the analytical results (Fig. [1} . Hence 
our numerical approach allows for a precise detection and 
location of any putative AT line. 

Results in three dimensions In the finite- 

dimensional system the scaling of the susceptibility re- 
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FIG. 3: (Color online) Finite-size correlation length as a func- 
tion of T/Tc for the 3D Ising spin glass at zero field with bi- 
modal interactions. The data cleanly cross at T c = 0.663(6). 

quires an additional parameter 77 whose putative in-field 
value is unknown, as opposed to the mean-field case. We 
thus prefer to study the transition via the scaling of the 
finite-size two-point correlation length 2(], 21 1 given by 



1 



2sin(|k min |/2) 



X(0) 



-1 1/2 



X(kmin) 



(2) 



where k m j n = (27r/L,0, 0) is the smallest nonzero wave 
vector and x(k) the wave- vector-dependent spin-glass 
susceptibility defined as 



3 ik-(Ri— Rj) 



(3) 



Note that we used a slightly different definition of the 
susceptibility than in Ref. [lfj. However, the two ex- 
pressions are equivalent, except that the definition used 
here re quir es only two replicas (unlike the four needed 
in Ref. [lfj) thus saving computational resources. The 
correlation length divided by the system size L has the 
scaling relation £, L /L = X (^/"[T - T C (H)]), where v is 
the usual critical exponent. When T = T C (H) data for 
different L cross signaling the existence of a transition. 

Data for the 3D model in zero field are shown in Fig. [3] 
A clear crossing point at T c = 0.663(6) is observed, thus 
signaling the presence of a spin-glass state for system 
sizes up to L = 16 and in agreement with previous re- 
sults [l9(. In contrast, data along the finite-field branch 
[h 3B (T) = 0.7(T - 0.3)] (see Fig. [TJ show no sign of a 
transition (see Fig. [4]). The inset shows a zoom of the 
data for small fields where the crossing points (arrows) 
between the different lines wander towards zero for in- 
creasing L. This shows that at temperature T = 0.3, the 
instability is below h < 0.006. This should be compared 
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FIG. 4: (Color online) Finite-size correlation length as a func- 
tion of hw(T) = 0.7 (T - 0.3) (finite-field branch in Fig. [TJ 
for the 3D Ising spin glass. The apparent crossing of the data 
wanders towards zero for increasing system size (see inset for a 
zoom into the low-field region), i.e., there is no sign of a tran- 
sition for H > (arrows mark crossings between L-pairs). 

with the mean field case where, e.g., for c = 3 with Gaus- 
sian couplings (having a similar T c as the 3D model) , one 
observes h(0.3) « 0.26, i.e., a value 40 times larger. Per- 
forming a one-parameter scaling with fixed volume ra- 



tio h c {L) ~ hf + b/L 2 gives h c c 



-0.00004(14), 



i.e. 



compatible with zero with a Q-factor of 94.9%. Fixing 
= gives Q = 99.7%; i.e., the data suggest the ab- 
sence of a spin-glass state in a field. 

Summary and discussion - We have presented cal- 
culations of the AT line of a spin glass on a Bethe lattice 
using the cavity method. These results for the mean-field 
model on the Bcthe lattice allowed the validation of our 
Monte Carlo simulations where we have observed an AT 
line close to its predicted value. This is in stark contrast 
to the 3D Ising spin glass with bimodally-distributed 
bonds where data for a wide range of fields and tem- 
peratures clearly show a lack of ordering in a field. 

It would be interesting to study higher-dimensional 
systems using the method presented here to verify 
whether or not an AT line is observed above the upper 
critical dimension 12j and we hope our results will also 
spark new theoretical developments in this direction. 
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